Dynamics of test bodies with spin in de Sitter spacetime 
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We study the motion of spinning test bodies in the de Sitter spacetime of constant positive 
curvature. With the help of the 10 Killing vectors, we derive the 4-momentum and the tensor of spin 
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one needs to impose the so-called supplementary condition. We discuss the dynamics of spinning 
test bodies for the cases of the Frenkel and Tulczyjew conditions. 
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I. INTRODUCTION 

The study of motion of test bodies in General Relativ- 
ity (GR) theory has a long history, see for example the 
account in Among other approaches, the multipole 
approximation method represents a powerful technique 
with the help of which one can derive a self-consistent 
set of equations of motion for a body characterized by 
the moments of arbitrary order. In the zeroth - or pole 

- order one recovers the geodesic equation in the context 
of multipolar methods. At the next - first or pole-dipole 

- order, the test body is described by the 4-momentum 
and the tensor of spin, and the dynamics is governed by 
the Mathisson-Papapetrou equations, the relevant dis- 
cussion can be found in }2-l20l|. Even in the absence of 
spin the integration of the equations of motion is a diffi- 
cult problem for nontrivial spacetimes. The spinless test 
body moves along a timelike geodesic on the curved man- 
ifold. When the spin is non-zero, the motion becomes 
nontrivial even in the fiat spacetime [HI, [Hj]. Moreover, 
in curved spacetime the motion is no longer geodesic due 
to the Lorentz-like force that acts on the test body. The 
Mathisson-Papapetrou force depends on the spacetime 
curvature and this considerabl y c omplicates the problem 
of finding the trajectories, [23H28l |. 

Here we investigate the dynamics of test bodies with 
spin in de Sitter spacetime. The latter is the maximally 
symmetric 4-dimensional space which means that there 
exist 10 (= 4x(4+l)/2) Killing vector fields that describe 
the symmetries of this manifold. The flat Minkowski 
spacetime also has 10 Killing vectors, and in this respect 
the geometrical properties of the de Sitter spacetime are 
close to those of the Minkowski space. At the same time, 
the de Sitter manifold has a nontrivial curvature. 

We should clearly stress at this point, that this work 
is not concerned with the actual derivation of the equa- 
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tions of motion in the context of different multipolar ap- 
proximation schemes. We will not discuss the conceptual 
questions which eventually lead to the imposition of dif- 
ferent supplementary conditions, nor do we discuss dif- 
ferent "flavors" of multipolar schemes - in particular not 
the subtleties regarding the definition of the moments 
within these schemes. Here we are only concerned with 
the solution of the pole-dipole equations motions in a spe- 
cific background spacetime for two frequently used sup- 
plementary conditions 0, H, HJ. As our analysis will 
show, different choices of the supplementary condition 
lead to quite different dynamics. In other words, the se- 
lection of such a condition - which in the context of the 
multipolar schemes under consideration has the status 
of additional assumption on the level of the equations of 
motion - should be performed with utmost care. These 
assumptions should not be misunderstood as the initial 
conditions. The number of equations of motion is less 
than the number of dynamical variables (see the next sec- 
tion) , and the supplementary condition makes the system 
self-consistent and predictive. 

We should also stress that our analysis is valid for both 
interpretations of the Mathisson-Papapetrou equations 
which can be found in the literature, i.e. it applies to 
point particles [29H32j as well as to extended test bod- 
ies. Recall that - on the level of the equations motion 
- the descriptions of both types of objects formally co- 
incide. One should keep in mind though, that there are 
preferences regarding the supplementary condition (33r - 
|35[, depending on the system which is supposed to be 
described by the equations of motion. 

The structure of the paper is as follows. In section |TT] 
we present a short overview of the equations of motion 
for spinning test bodies. We pay particular attention to 
the conserved quantities in these equations. In section 
Mil we collect some facts about the Killing vectors of dc 
Sitter spacetime. These results are then used for the 
integration of the equations of motion in de Sitter space- 
time in section IIV1 This is followed by a discussion of 
how momentum and spin can be expressed with the help 
of the integrals motion in section |Vj Finally, we draw 
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our conclusion in I VII and present a brief outlook on open 
problems. A summary of our conventions and a directory 
of symbols can be found in the appendix [A] 

II. EQUATIONS OF MOTION 

A test body with spin is described by the following 
variables. Its trajectory is given by 4 spacetime coordi- 
nates x a (s) as functions of the afhne parameter s (proper 
time). Furthermore, the body is characterized by the first 
two moments (pole and dipole) coming from its energy- 
momentum contents: the 4-momentum p a and the spin 
S a P = — S^ a . One can interpret these as the two "grav- 
itational charges" carried by the body since they couple 
to the gravitational held in the same way as the electric 
charge couples to the electromagnetic held. 

The polc-dipolc equations of motion (usually known as 
the Mathisson-Papapetrou equations) read 

P a = -\s^u?R^ a , (2.1) 
fr*P _ 2p[ Q iA (2.2) 

Here u a = dx a /ds is the 4- velocity of the body, and the 
dot denotes the covariant derivative with respect to the 
proper time, " ' " = D/ds = u a V a . The force term 
on the right-hand side of (12.11) depends explicitly on the 
spacetime curvature. 

The set of the Mathisson-Papapetrou equations is in- 
sufficient to determine the dynamics of the system. In- 
deed, we have 14 unknown variables (x a ,p a ,S liU ) and 
only 10 equations (|2.ip - (|2.2p . One additional algebraic 
equation comes in the form of the normalization condi- 
tion u a u a = 1. Thus one needs 3 more equations to 
make the model predictable. Such additional equations 
are usually algebraic and they are commonly known as 
the supplementary conditions. Although this name is 
rather misleading, we will keep the tradition. The two 
most widely used are the Frenkel condition [2l[ (some- 
times also called Pirani condition Q) 

S°% = 0, (*) (2.3) 

and the Tulczyjew condition Q 

S a Pp = 0. (**) (2.4) 

Both algebraic equations have 3 independent components 
and thus the total number of the equations becomes equal 
to the number of the unknowns. The Frenkel condition 
(I2.3[) was introduced in the electrodynamical context by 



Since the dynamics of a spinning particle is non-geodesic, in gen- 
eral, this cannot be viewed as an integral of motion, but is a mere 
constraint due to the choice of paramctrization of the trajecto- 
ries. 



taking into account the purely "magnetic-dipole" nature 
of electron's spin, and it is nowadays mainly used in a 
point particle context. Although also Tulczyjew's condi- 
tion (12.41) was initially used in a point-particle representa- 
tion, it is nowadays more often associated with extended 
bodies. 

In this paper we will analyze the dynamics of spinning 
test bodies for both supplementary conditions. 
By contracting (|2.2[) with up, we hnd 

p a = mu a + S afi Uf}. (2.5) 

The m := u a p a we will call the rest mass of the body, 
defined as usual as the projection of the 4-momentum on 
the rest frame of moving body. Besides that, we can de- 
fine another mass parameter by M 2 := p a p a - In general, 
these two masses arc different. We will compare them be- 
low. Depending on the supplementary condition chosen, 
the mass parameters may be constant or not. 

A. Conserved quantity 

Let £ be a Killing vector. This is a solution of the 
equation V a £/3 + V^£ a = 0. Applying V 7 , we derive 

V 7 V Q £/3+V 7 V^a = V 7 V Q e/3+V /3 V 7 Ca-i? 7 ^ A 6 = 0. 

Now, add and subtract V Q V/j£ 7 and use the identity 
V 7 V Q £,3 + V^V 7 ^ a + V Q V,3£ 7 = 0. As a result, we ob- 
tain the second covariant derivative of any Killing vector 
in terms of the curvature: 

V a Vp^ = Rp- (a x £x. (2.7) 

Then we straightforwardly hnd 

^(2£ aP Q ) = 2i a p a + 2^ a p a , (2.8) 

^(S Q/5 V Q &) = -2& s p a + S afl u*R a f h x £x- (2.9) 

In the last equation we used (12. 2|) and ()2.7|) . Taking the 
sum of (|2.8|) and (|2.9[) , and using the equation of motion 
(|2.1[) . we derive 

£ (2£ aP a + S a0 W a ^) = 0. (2.10) 
Thus, we have demonstrated that the scalar 

2£ a p a + S a(3 V a & = const (2.11) 

is conserved. Quite remarkably, no supplementary condi- 
tion is needed. For other conserved quantities, non-linear 
in spin, see [H, H3| • 

B. Mass parameters M and m 

Here we study whether the mass parameters are con- 
stant or may change along body's trajectory. Recall the 
definitions m — p a u a and M 2 = p a p a . 
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Equation (|2.1[) yields u a p a = 0. Consequently, 

m=p a u Q = u a £-(S afi up). (2.12) 

In the last equality we used (|2.5|) . As we see, m is con- 
stant when the Frenkel condition (|2.3p is assumed. 
If we contract (|2.2p with p a p,g, we find 

pp — (p a S al3 ) + mMM = 0. (2.13) 
Accordingly, M is constant for the Tulczyjew condition 



C. Velocity and momentum 

The relation between the velocity of the body and its 
momentum provides a very useful information for the 
integration of the equations of motion. Here we derive 
this relation explicitly without using the supplementary 
conditions. 

Contracting (|2.2I) with pp and making use of (|2.1[) . we 
find 



1 D 



■P 



1 



M 2 * M 2 ds v ^ p > 2M 2 

(2-14) 

Now let us derive a useful identity. For an arbi- 
trary skew-symmetric tensor (p a p we define the dual 
._ l^aPftu p^ v . Any two such tensors together with 
their duals satisfy the identity (a somewhat lengthy but 
straightforward proof is based on the properties of the 
Kronecker and Levi-Civita objects, see [38| chapter 7, 
e.g.) 



p ai ap 1 



r Vp 1 = 2 °P v a > 



ft is • 



Taking a — (p* and noticing that if* 
from (I2.15P a new identity 



f 1 V/3 1 = -7 Op Iff tp^. 



(2.15) 
(p, we obtain 

(2.16) 



These identities underlie the following important alge- 
braic result. Let us consider the second rank tensor 



K a p :=,5f + ^r /37 . 



(2.17) 



Here p ai and ap 1 are the two arbitrary skew-symmetric 
tensors. Then the inverse of (j2.17[) reads 



(K- 1 ) a p = - 



[i + |(^)] 6<$-<fn<7fr 



(2.18) 



+ i(^)-n(w*)( f7(7 *) 

Here ((pa) = p^a^ and ((pip*) = p^p*^. The proof is 
straightforward: one should multiply (K~ v ) a 1 K~ < p and 
make use of both identities (|2.15p and (|2.16[) . 



We can write f|2 . 14[) as the algebraic system 

Ka ^ = > a ~w^ sa ^ (2 - 19) 

where p afi = S af} and a a p = ^S^R pua p. Accord- 
ingly, we find (pa) — ^^S^ S a ^ R^ap, for example. 

Then we find that the velocity of a test body can be 
expressed in terms of its momentum and spin: 



u a = (K~ 



n0 



771 „«* 1 D {Qctfi n \ 



P p > p 

(2.20) 

This result is valid for any supplementary condition. For 
the case of the Tulczyjew condition this relation was 
originally derived in |39l - l41| . When spin satisfies (any 
of) the supplementary conditions (|2.3[) or (|2.4j) . we have 



((pip*) = S^S, 



a/3 



0, and hence 



2 g a f> SP ,v R 



[ivfS~fP 



AM 2 + S a ?S^R 



(2.21) 



III. DE SITTER SPACE: KILLING VECTORS 

Let -da be a coframe 1-form. The curvature of the de 
Sitter spacetime reads 



Rj = ^ a A^. 



(3.1) 



Here £ is a real constant. The anti-de Sitter space arises 
with the help of the formal replacement I 2 ~ > — l 2 . In 
components, R a ^ = ^R^a^^ A we have R^a 13 — 
ji (g a ^ - g av S^). 

The de Sitter spacetime has many faces. Depending on 
the choice of the local coordinates, the metric can have 
static form, either isotropic or Schwarzschild-like, or it 
can be written in the cosmological form of an expanding 
world. 

In static isotropic coordinates, the line-element of the 
de Sitter spacetime reads 

ds 2 = V 2 dt 2 - W 2 (dx 2 + dy 2 + dz 2 ). (3.2) 

The functions depend only on r 2 = x 2 + y 2 + z 2 : 

l-r 2 /£ 2 „. 2 



V 



l + r 2 /£ 2 



W 



l + r 2 /( 2 



(3.3) 



The coordinate transformation X a = Wx a , a = 1, 2, 3, 
(x 1 = x, x 2 = y, x 3 = z) hence 



p = W r = 



2r 



l + r 2 /£ 2 ' 



(3.4) 



brings the line element to the standard spherically sym- 
metric form 



ds J = 1 - hr }dt 



dp 2 



- p 2 dd 2 -p 2 sin 2 9 d<p 2 



1 p 
l ^ -p 



(3.5) 
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Here p 2 = X 2 + Y 2 + Z 2 and spherical coordinates 
are introduced in the usual way by X a — {X — 
psinflcos^, Y — p cos 8 cos (f>, Z = pcosd}. The met- 
ric (|3.5J arises from the Kottler (Schwarzschild-de Sit- 
ter) |42l . metric when the mass of the central source 
is zero. 

Another change of coordinates from the static to the 
cosmological frame 



t 

X a 



t + -\ og (i-p 2 /e 2 ) : 



X a 



V 1 - p 2 /p 



(3.6) 

a = 1,2, 3, (3.7) 



brings the de Sitter metric into the form of an exponen- 
tially expanding world 



ds 2 = dt 2 



S ab dX a dX b 



(3.8) 



The de Sitter spacetime can be viewed as a hyper- 
boloid embedded into the fiat 5-dimensional spacetime. 
Let X A , A = 0, 1, 2, 3, 4, be the coordinates and the line 
element ds 2 = r)ABdX A dX B of such a spacetime with 
t)ab = diag(l, — 1, — 1, — 1, — 1). The de Sitter manifold 
can then be embedded in it as the hyperboloid 

(X ) 2 - (X 1 ) 2 - (X 2 ) 2 - (X 3 ) 2 - (X 4 ) 2 = -I 2 , (3.9) 

where the embedding coordinates can be chosen for ex- 
ample as 



l-^_sinh(</£), 



X y] 



• V = >' i / I ~ ~p cosh(V4 
X a = X a , a= 1,2,3. 



(3.10) 

(3.11) 
(3.12) 



A. Killing vectors 

The de Sitter manifold is a maximally symmetric 
space, and there are 10 Killing vector fields on it. In 
the isotropic coordinates {t,x a ) 1 the Killing vectors read 
explicitly: 



r(0) 



fit, 



£(«) = e abc x b — , 



d_ 

' dx c 



(a) 



±t/l 



W x a d 
VTdi 



± 



V 



•ab 



e 2 



d 
dx 1 



(3.13) 
(3.14) 

.(3.15) 



In the static coordinates (t,X a ), they look very similar 

fa 



£(0) 



^abc v 

e A b 



d 



(a) 

v± = 



dX cl 
e ±t/t 

V 1 - p 2 I (2 



Zl±± (i -P- 

e m \ p 



8X a 



(3.16) 
(3.17) 

.(3.18) 



In the cosmological setting, the Killing vectors read 

(3.19) 
(3.20) 

V+ : — • (3-21) 

V- 



dt I 8X a ' 
= e abc X b - 8 



dX c 



» 



d 



dX a 

(a) _ 2X a f d X» d 



I \dt i dX b 



P_ _ -It /I 

I 2 



d 



dX a 



Here p 2 = ^X a ) . The expression for r^L 1 is rather non- 
trivial. 



(3.22) 



B. Conformally flat representation 

Since the Weyl tensor for the de Sitter space is trivial, 
the metric can be recast into a form that is conformally 
flat. Explicitly, 

ds 2 = tp 2 r\ijdx l dx\ rjij — diag(+l, — 1, —1, —1). 

(3.23) 

The conformal factor depends only on the 4-dimensional 
"radius" a — rjijX l x'J , namely, 



1 



ip = 



1 



ir 2 



1 



IP 



(3.24) 



In this representation, the Killing vectors are as follows: 



1+ 4£ 2 ) da U 2 



(a) 

£ = x a dp - xpd, 

[a/3] 



dp, (3.25) 
(3.26) 



IV. INTEGRATING THE EQUATIONS OF 
MOTION IN DE SITTER SPACETIME 



In de Sitter spacetime with the curvature (|3.ip . the 
equation of motion (|2.ip reduces to 



p ~ g2 a U P- 



(4.1) 



The complete integration of the dynamical equations de- 
pends crucially on the supplementary condition. The two 
most important cases are analyzed separately below. 



A. Tulczyjew condition 

Assuming (|2.4[) . we introduce the 4- vector of spin via 
S a := r\ a t ^ l 'p >pS 'nu ■ The inverse formula yields the spin 



5 



tensor in terms of the spin vector: S a ^ = j/gr^"' 3 '" 'p^S v . 
By construction, we have the orthogonality properties 



Pa S a = 0, S aP S a = 0. 



(4.2) 



With the help of (|4.1I) and (|2.14p we derive further or- 
thogonality properties 



i>J a = o, 



i a S n = o. 



(4.3) 



An immediate consequence is the (covariant) constancy 
of the spin vector. Indeed, we find 



S a = 



1 



2NP 



p a ppSP)=0. (4.4) 



Using the orthogonality properties, we find in the de Sit- 
ter space 

= ^2 (P a P?S 2 ~ 6 P M2 S 2 + M 2 S a S p ) .(4.5) 
Consequently, 

S 2 



2M 2 P \M 2 
Substituting this into (I2.14p . we have 

S 2 



1 - 



4M 4 £ 2 I \M 2 



(- 



P 



= 0. 



(4.6) 



(4.7) 



Thus, since the spin vector is spacelike, we can never have 
S 2 = 4M 4 £ 2 , and hence we obtain 



P 



m = M. 



(4.8) 



As a result, the equations of motion of a spinning test 
body in the de Sitter spacetime under the Tulczyjew con- 
dition reduce to 



r ( => o, 

or, equivalcntly 



a (**) o 

p = mu 



0. 



o. 



(4.9) 



(4.10) 



The first equation actually means that the trajectories 
of the spinning bodies are the geodesies in the de Sitter 
space. The second equation describes the precession of 
the spin vector, or tensor, of a body during its motion 
along a geodesic curve. 



Frenkel-Pirani condition 



have a certain similarity to the above case, there are im- 
portant differences. In particular, from (|4.1|) it immedi- 
ately follows that, like in the previous case, the momen- 
tum is covariantly constant, p a = 0. 

Following the same line of reasoning, we define the 4- 
vector of spin by S a := rfP^upS^. The inverse for- 
mula yields the spin tensor in terms of the spin vector: 
gap _ i^a^nuy^^g^ j- we uge ^ ne norma li za tion u 2 = 1). 

By construction, we thus have the orthogonality proper- 
ties 



u a S a = 0, S af3 S a = 0. 



(4.11) 



Now in complete analogy with (|4.4[) . directly from the 
definition of the spin vector, we derive that the spin vec- 
tor is Fermi- Walker transported: 

S a = \v a ^ v u^ vpi7 uPS a = -u a u^. (4.12) 



We thus have the system 2 

P a = 0, p^S f3 = 0. 



(4.13) 



Although this looks formally similar to (|4.9p . the actual 
dynamics is very different. In particular, the trajectories 
are no longer geodesies because the momentum does not 
coincide with the velocity. Instead, 



P 



S aP v 



up, 



(4.14) 



and this relation must accompany the integration of the 
system (|4.13j) . The above system can be simplified even 
further. If we differentiate eq. (14.141) and contract it with 
S a , we obtain - with the help of (|4. 1 1|) : p a S a — mu a S a — 
S a PupS a . The last term vanishes when we use (|2.2p . 
and the left-hand side vanishes because of the covariant 
constancy of the momentum, cf. (|4.13p . Thus, we obtain 
u a S a = 0, in other words - taking into account (14. 12|) - 
the spin is also parallely transported in the Frenkel-Pirani 
case. We thus end up with the final system: 

p a ( = } 0, S a = 0, p a ( = } mu a - S at3 Uf). (4.15) 
Equivalently, one may look for solutions of the system 

S al3 up-mu a ( = } 0, (4.16) 

(4.17) 



S al3 + 2u [a S^ii 1 = 0. 



The geodetic motion plus the parallel transport of the 
spin 



0, 



S a = 0, 



(4.18) 



is a solution of (|4.15l) . as one can immediately check. 
However, in general a spinning body does not move along 
a geodesic. 



Let us now analyze the Frenkel case (12. 3p . Although 
the dynamic equations for this supplementary condition 2 Here 
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The force that pushes the body away from a geodesic 
is produced by its own spin, and the resulting motion 
is a classical analog of the Zitterbewegung. The key to 
the description of this motion is a new vector variable 
that can be introduced along the lines of the fiat space 
discussion 14511 . i.e. 



(4.19) 



In view of the covariant constancy of the momentum, 
both invariants M 2 = p a p a and 



(4.20) 



are conserved along body's trajectory. As a result, the 
length of the new vector field is also constant 



(4.21) 



Being orthogonal to the velocity, Q a is spacelike, and 
thus it is rotating thereby generating the helical motion 
of the body. This can be further clarified as follows. 
Contracting (I2.2p with pp, we find 



mp 
~M 2 



0. 



(4.22) 



By differentiating we thus prove that the acceleration is 
produced by the "Q-force" 



(4.23) 



On the other hand, every vector can be expanded with 
respect to the natural orthogonal basis formed by the 
quadruple (p a , S a , Q a , Q a ) . This basis is not orthonor- 
mal since the lengths of the legs are not equal ±1, but one 
can straightforwardly verify that each vector is orthogo- 
nal to the three others. In particular, we can expand the 
acceleration with respect to this basis: 



ap a + (3S a + 1 Q a +8Q a . 



(4.24) 



Contracting this with p a ,S a , Q a , we find a = j3 = 5 = 
since u a p a = ii a S a — ii a Q a = 0. The only nontrivial 
coefficient is thus 



u a Qc 

Q 2 



(4.25) 



However, contracting (|4.14l) with p a , we find u a Q a = 1 — 
m 2 /M 2 , and then using (|4.21J) . we have 7 = - 4M 2 /S 2 . 
Accordingly, (|4.24|) reduces to 



4M 2 



(4.26) 



Comparing (|4.23[) with (|4.26p . we derive the oscillator 
equation 



Q a + co 2 Q a 



0. 



(4.27) 



The frequency is defined by 
2M 



M 



V^S^ R 



\lhs aP s^' 



(4.28) 



Furthermore, substituting u a from (|4.22|) into the equa- 
tions of motion and the Frenkel condition, we can recast 
(EH) and (EH) into 



a/3 = o 



^Qp-mQ a = 0. 
Here we introduced another interesting object 



(4.29) 
(4.30) 



;= Sa +paQ P_^ Qa = Sa 0_PP ls y0_P^ Lsali 

(4.31) 

that is the projection of spin on the momentum. Using 
it, the frequency (|4.28p is recast into 



rn 



(4.32) 



Qualitatively, the dynamics of spinning bodies subject 
to the Frenkel condition in the de Sitter spacetime is sim- 
ilar to that in flat space [45j]. Everything is determined 
by the initial conditions. If initially (at the proper time 
s = 0) spin is parallel to the momentum, i.e. S al3 pp = 
(hence Q a = 0), then this is true on the whole trajectory 
that turns out to be geodesic. Otherwise, the trajectory 
is a geodesic curve, perturbed by the oscillatory motion 
of Q a with the frequency (j4~28| and (|4T32|) . 



V. USING THE INTEGRALS OF MOTION 

As a matter of fact, the de Sitter spacetime has exactly 
the same number of Killing vectors as the total number 
of the "gravitational charges" , that is, 10. Then we can 
try to find the momentum p M and the spin S^ 11 without 
solving differential equations by just making use of the 
10 conservation laws. 

This task is most straightforwardly treated in the con- 
formally flat representation. By substituting (|3. 251) and 
(|3.26[) into (|2.11[) . we have the algebraic system 

+ = 2n Q , (5.1) 

(a) (a) 

2 e^ + ^V^ £„ - 2£ a/J . (5.2) 

[a/3] [a/3] 

Here n a and E Q( 3 = — are the 4 + 6 = 10 constants 
of motion. We introduced the factors 2 for convenience. 
It is worthwhile to notice that due to the skew symmetry 
of spin, S^Vptv = S^d^. 

For the translational (|3.25p and rotational (|3.26p 
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Killing vectors we find 



With the help of (pTT5|) and (f5TT6|) . we recast (|5T3)) into 



(a) (a) 
(a) 

[a/3] [a/3] 



2<p 4 

12 



(7 \ Xj/iCa 



"*^[/i^?^]a i 



[a/3] 

Here we introduced 

'/mo •- yi 4 ^ 2 j ^ + • 

In addition, we define a similar object 
_ /-. _o_\ 

»7^« •- + 4 ^ 2 ^ 2 ^ 2 • 



,(5.3) 

(5.4) 
(5.5) 
(5.6) 

(5.7) 
(5.8) 



It is easy to check that up to a factor they are inverse to 
each other, 



( 1 -4^)( 1 + ^)^ ™ 
Other useful relations are as follows: 



1 + _E_y 5 *-^ 

API v P 



(5.10) 
(5.11) 



The indices are everywhere raised and lowered with the 
help of the flat Minkowski metric t]^. 

Substituting (|B3] ) -([5T6 | into (f5Tj) and (|53|) . we derive 

4 

V 2 Va^P^ + -p f) av x ^S 11 " = II Q , (5.12) 

V 2 {ZaVPn - XpVa^P^ + <P 4 flafiVPv S^" = E a p. (5.13) 

As a first step, contracting f|5 . 1 2[) with r) al3 and using 
(EU) and (jSTTU)) . we find 



As a result, we have 

V 2 (XaVPM ~ x PVan)P >> 



On the other hand. 



= V 4 f) a nr}p»S flv . (5.16) 



S^Xu + fT^u). (5.14) 



-ir 2 



Xallp - Xpll a 

(5.15) 



" J ^CW" ^ = ~ 1 - a (^aII/3 - X^IIq) 



1 + 32* 



1 



(5.17) 

Now contracting with fj ap fj°' 7 and making use of (|5.9p . 
we find the spin tensor explicitly 

5"" = ^T" Sa/J + ^ Q n Q ic" - r) va n Q (5.18) 



Notice that 



— 1 



XT, 



(5.19) 



which we took into account when deriving (|5.18l) . Using 
this again, we find 



S^Xu = 



i r 



(l + ^) rT^apxP + W v o - x»x») U> 



(5.20) 

Substituting this into (|5.14[) we finally obtain the mo- 
mentum 



1 



P' 



(5.21) 



In summary, by using the 10 first integrals correspond- 
ing to the Killing vectors of the de Sitter spacetime, we 
are able to express the momentum and the spin as func- 
tions of the constants of motion: 



P — ^2 Z-'afiX -f A/ 



(5.22) 



5^ = fj fia fj^T, a p + ^ a Il a x v -r) va Il a x^.{5.23) 

Remarkably, the dependence on the spacetime coordi- 
nates is merely polynomial. Notice that two different 
etas appear in the final formulas: ff v enters (|5.22p but 
rf v enters ([5~35]) . 

This seems to be the best what one can do without im- 
posing the supplementary condition on spin. The dynam- 
ics of momentum and spin is completely known, given by 
(|5.22[) and (|5.23[) . however, the trajectory of the body is 
still undefined. At first sight, one may think that it is 
possible to substitute spin and momentum into (|2.5|) and 
solve the resulting algebraic equation for the 4-velocity 
u a . When this is done, one can find the trajectory from 
the velocity vector field. But this plan does not work be- 
cause one can verify by direct substitution of (|5.22[) and 
(|5.23[) into (|2.5[) that the latter is an identity. 

After imposing the supplementary condition, every- 
thing reduces to a mere technical problem of integrat- 
ing the first order system. For the Tulczyjew condition, 
such a system is obtained by substituting (|5.22[) into the 
left-hand side of (|4.8[) and recalling that u a = dx a /ds on 
the right-hand side. Similarly, for the Frenkel condition 
one needs to plug (|5.22[) into (|4.22[) . with an intermediate 
step of constructing Q a by contracting (|5.22p with (|5.23p 
and substituting the result into (|4.22p . In both cases, for 
the Frenkel and the Tulczyjew condition, the resulting 
first order system involves only polynomial functions of 
x a and can be straightforwardly integrated numerically. 



8 



VI. DISCUSSION 

We have studied the dynamics of spinning test bod- 
ies in the de Sitter spacetime. Qualitatively, the re- 
sults are similar to those obtained in flat spacetime. For 
the Tulczyjew condition (12.41) . the body moves along the 
geodesic curve, whereas the spin vector is parallelly trans- 
ported along the trajectory. In the Frenkel case (|2.3I) . the 
spin is still parallelly transported, but geodesic motion 
is just one special solution of the equations of motion. 
When the initial value of Q a = S^pn/p 2 is nontrivial, 
then the body is affected by the spin-dependent force, the 
acceleration u a is nontrivial, and the trajectory oscillates 
around a geodesic with the frequency (14.281) . (I4.32[) . The 
curvature of the de Sitter space thus affects the dynamics 
only indirectly through the structure of the correspond- 
ing geodesies on this manifold. The introduction of the 
variables fi a p and Q a can be qualitatively compared to 
the definition of the mean spin and the mean position 
operators in quantum mechanics (46l - |49| . 

In summary, different supplementary conditions do 
lead to fundamental changes on the level of the equations 
of motion. In particular we have explicitly shown for the 
de Sitter spacetime, how the solutions for the worldlinc 
change under the Tulczyjew (|2.4p and the Frenkel (|2 .3[) 
condition. The search for further solutions of the mul- 
tipolar equations of motions is an ongoing task, in par- 
ticular one should aim for a solution of the equations in 
more complicated spacetimes than the one covered in the 
present work. 

As mentioned in the introduction there also remain 
conceptual questions to be answered in the context of 
different multipolar approximation schemes. In particu- 
lar the interpretation of the quantities and the structure 
of the final set of equations of motion in such schemes 
should be investigated. Without going into detail at this 
point - see the upcoming work [50( in this respect - we 
have to stress that the system (|2.1I) - (|2.2I) can be obtained 
in several different ways. While one can achieve formal 
equivalence on the level of the equations of motion, in 
particular up to the pole-dipole order, there are differ- 
ences in the derivation as well as in the interpretation of 
these equations. 

Let us briefly mention two viewpoints here. Most in- 
terestingly, the equations (|2.1[) and (I2.2p have been in- 
terpreted in a point-particle, as well as in an extended 
body context. In the point-particle picture it is imme- 
diately clear that the dynamics of the particles under 
consideration is directly influenced by the choice of the 
supplementary condition. The particle is thought to be 
localized at the worldline, and the worldline represents an 
immediate description of its motion through spacetime. 
This should eventually lead to small, but physically de- 
tectable changes in the motion of these particles. For 
an indirect detection one may imagine a charged particle 
moving in an electromagnetic field. The differences in its 
motion due to the different supplementary conditions - 
recall the oscillatory motion in case of the Frenkel con- 



dition - should be detectable via radiation losses. For 
an extended body the impact of different supplementary 
conditions may be less straightforward to detect. The 
supplementary condition can be interpreted as the choice 
of a suitable representative worldline of the body under 
consideration - in analogy to the dynamics of bodies in 
Newtonian physics one may think of an extension of the 
concept of the center-of-mass. While it is clear that the 
shape of this worldline will change under different sup- 
plementary conditions, the question remains up to which 
level such a change impacts observable quantities. Micro- 
scopic oscillations along the representative line probably 
do not play any role in the description of the motion of, 
say, an extended star. Nevertheless one has to check for 
each application, if the length scales of the system under 
consideration justify to view the choice of the supplemen- 
tary condition as something which has no direct physical 
impact. 

Let us close by pointing out, that further work on the 
foundations and the interpretation of multipolar schemes 
is needed. This concerns both, i.e. the point-particle as 
well as extended body, interpretations. In particular, one 
should always keep in mind that the multipolar schemes, 
which lead to equations of motion (|2.ip and (12.2[) . are 
approximation schemes which - by construction - only 
capture certain features of the full theory. 

Although the present study is mainly of theoretical 
nature, one should keep in mind that a deeper under- 
standing of the dynamics of spinning bodies in curved 
spacetime is important for many physical applications. 
In particular, the spin effects are well observable in as- 
trophysical situations for the binary star systems [5~iT - [53j . 
Furthermore, the use of the conserved quantities (I2.11[) 
appears to be quite useful also in situations when the 
number of Killing vectors is insufficient for the complete 
integration of the equations of motion, cf. for example 
the study of spinning bodies in the gravitational fields of 
the black body type sources in [55| . 
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Appendix A: Dimensions & Symbols 

In order to fix our notation, we provide some tables 
with definitions in this appendix. The dimensions of the 
different quantities appearing throughout the work are 
displayed in table |U Table [XT] contains a list with the 
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TABLE I: Dimensions of the quantities. 



Dimension (SI) 



Symbol 



most important symbols used throughout the text. Greek 
indices denote 4-dimensional indices and run from a — 
0, . . . , 3, the signature is (+,-,-,-). 



Geometrical quantities 



(«) 

x - dx - s , 



1 a p ' 

Rap^f , Ra/z 



Matter quantities 



1 

kgm/s 
kg m 2 /s 
kg 2 m 3 /s 2 
kg/m 2 s 



m, M, p a , H a 

Qab no lt ab v 

s a 

rpaP 



Auxiliary quantities 



m 

™2 



K a p, W, V, (fi, f) a p, fj a 

£, p, Q" 

a 

u> 

7 



Operators 



1 


Pp-i * : V 


m" 1 


r) V £.—«■» 

Wl ' V " ds - 
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_D _«■ v 
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